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SUMMARY 


A  method  for  computing  the  geometry  of  the  tip  vortex  of 
a  helicopter  rotor  in  steady,  forward  flight,  including  the  dis¬ 
tortion  due  to  the  velocities  induced  by  the  vortex  wake,  has 
been  developed.  The  relative  importance  of  various  parameters 
and  approximations  used  by  this  method  has  been  evaluated  for  one 
test  case,  including  comparison  with  experimental  results  ob¬ 
tained  from  a  flow  visualization  study. 
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I .  INTRODUCTION 


As  rotary  winged  aircraft  find  wider  applications  greatly 
improved  performance  is  being  required  of  them,  e.g.,  greater 
speed,  better  aerodynamic  efficiency,  and  lower  vibration  levels. 
To  meet  these  requirements,  a  better  understanding  of  rotor  aero¬ 
dynamics  is  needed. 

A  recent  major  improvement  has  been  the  replacement  of  the 
classical  uniform  inflow  assumption  with  the  uniform  inflow  due  to 
the  tilt  of  the  tip  path  plane  (TPP)  relative  to  the  free  stream 
plus  the  nonuniform  induced  inflow  cue  to  the  vortex  wake  of  the 
rotor.  This  induced  inflow  due  to  the  vortex  wake  has  been  cal¬ 
culated,  with  the  aid  of  high  speed  digital  computers,  by  numeri¬ 
cally  integrating  the  Biot-Savart  relation  for  the  velocity  in¬ 
duced  by  a  vortex  element  over  some  assumed  vortex  wake  geometry 
(e.g.,  see  Refs.  1  through  4) .  The  most  common  vortex  wake  geo¬ 
metry  assumption  has  been  a  "rigid  wake".  This  assumes  that  the 
vortex  wake  has  the  spiral  shape  traced  out  by  a  point  on  one  of 
the  rotor  blades  parallel  to  the  tip  path  plane  (TPP)  and  drifts 
perpendicularly  to  the  TPP  with  the  classical  uniform  inflow  ve¬ 
locity.  Reference  1  also  tried  a  more  refined  approach  which 
left  the  rigid  wake  assumption  unchanged  parallel  to  the  TPP,  but 
calculated  the  nonuniform  induced  inflow  at  the  TPP  using  a 
rigid  wake  assumption  and  used  this,  or  some  of  its  harmonics, 
plus  the  uniform  inflow  component  of  the  free  stream  to  determine 
the  drift  velocity  of  the  vortex  wake  perpendicular  to  the  TPP, 

The  approach  taken  here  is  to  assume  a  rigid  wake  and  cal¬ 
culate  the  time  history  of  induced  velocity  at  various  points  in 
the  wake.  This  is  integrated  to  obtain  an  improved  estimate  of 
the  wake  geometry.  This,  in  turn,  is  used  to  calculate  a  new 
time  history  of  induced  velocity,  etc.,  and  the  iteration  con¬ 
tinues  until  convergence  is  achieved.  Various  refinements  to 
improve  convergence  are  required,  but  this  iterative  approach 
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still  requires  a  good  deal  less  computation  than  allowing  the 
wal'.e  to  distort  continuously  as  the  calculation  of  Induced  ve¬ 
locities  proceeds,  as  in  Ref.  9. 

This  report  starts  with  a  discussion  of  the  choice  of  a 
model  for  the  vortex  wake.  The  development  of  the  computational 
procedure  used  is  outlined  and  some  alternate  techniques,  now 
discarded,  are  discussed.  Some  preliminary  results  are  pre¬ 
sented  and  some  recommendations  are  made  for  future  development. 
The  appendix  presents  a  detailed  outline  of  the  computational 
procedure. 
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II.  VORTEX  WAKE  MODEL 


A  detailed  model  for  use  in  calculating  rotor  airloads 
would  represent  the  rotor  blades  by  lifting  surfaces  and  the  rotor 
wake  by  vortex  sheets.  The  calculation  of  the  wake  geometry  would 
involve  the  computation  of  the  distortion  and  roll-up  of  these 
vortex  sheets  due  to  their  own  induced  velocities  and  those  of  the 
lifting  surfaces.  As  the  vortex  sheets  rolled  up  into  line  vor¬ 
tices  viscous  effects  would  become  important  in  the  vortex  core. 

Unfortunately,  the  above  model  is  far  too  ambitious  for 
practical  computations  at  present  and  a  simplified  model  must  be 
constructed.  The  rotor  blades  can  be  represented  by  lifting  lines 
(bound  vortices)  instead  of  lifting  surfaces.  The  vortex  wake  can 
be  represented  by  a  set  of  trailing  and  shed  line  vortices  instead 
of  vortex  sheets.  The  trailing  vortex  lines  are  trailed  aft  along 
streamlines  relative  to  the  rotor  blade  while  the  shed  vortex 
lines  are  shed  parallel  to  the  instantaneous  position  of  the  rotor 
blade . 

For  maximum  thrust,  the  angle  of  attack  and  hence  the  bound 
circulation  should  be  as  nearly  constant,  radially,  as  possible 
over  the  outer  portion  of  the  blades,  where  most  of  the  lift  is 
generated.  An  idealized  model  of  the  radial  bound  circulation 
distribution  on  a  rotor  might  have  constant  bound  circulation 
from  the  50  percent  radius  to  the  tip  and  zero  circulation  in- 
ooard.  This  would  result  in  a  pair  of  line  vortices  being  trailed, 
one  from  the  tip  and  one  from  the  50  percent  radius. 

A  relatively  sharp  cutoff  of  the  bound  circulation  at  the 
tip  of  the  rotor  blade  is  normally  obtained  in  actual  rotor  opera¬ 
tion,  and  a  concentrated  tip  vortex  is  observed  (Ref.  5).  The  in¬ 
tersection,  or  near  intersection,  of  this  concentrated  tip  vortex 
with  a  rotor  blade  is  thought  to  be  a  major  cause  of  the  high 
harmonic,  almost  impulsive,  airloads  observed  in  flight  (Ref.  6). 
The  location  of  this  tip  vortex  is,  therefore,  likely  to  have  a 


considerable  effect  on  the  high  harmonic  airloads. 

The  inboard  bound  circulation,  however,  normally  tapers  off 
gradually  in  actual  rotor  operation,  thus  generating  a  broad 
vortex  sheet.  It  is  questionable  whether  this  inboard  vortex 
sheet  ever  rolls  up  tightly  to  become  as  concentrated  as  the  tip 
vortex.  In  any  case,  it  does  not  seem  to  roll  up  near  enough  to 
the  rotor  to  induce  large,  high  harmonic  airloads  similar  to 
those  induced  by  the  tip  vortex.  Since  the  determination  of  the 
distortion  and  roll-up  of  the  inboard  vortex  sheet  requires  too 
much  computation  for  practical  applications,  and  since  the  geome¬ 
try  of  this  inboard  vortex  sheet  does  not  affect  the  airloads 
nearly  as  much  as  that  of  the  tip  vortex,  the  inboard  vortex 
sheet  is  represented  by  a  set  of  undistorted,  rigid  wake  vortex 
lines. 

Experience  with  rotor  airload  calculations  has  shown  that 
for  the  far  wake  (the  part  of  the  wake  more  than  ^5  deg  to  90  deg 
behind  the  rotor  in  azimuth),  a  two  trailing  vortex  model  (a  tip 
vortex  plus  one  inboard  vortex)  gives  essentially  the  same  results 
as  a  many  trailing  vortex  model  since  the  circulation  is  sub¬ 
stantially  constant  along  the  blade.  For  the  near  wake  (the  first 
45  deg  to  90  deg  in  azimuth  behind  the  rotor  of  the  wake)  a  many 
trailing  vortex  model  is  often  necessary  in  airload  calculations 
when  a  rapid  change  in  circulation  exists.  However,  for  a  wake 
geometry  calculation  induced  velocities  are  calculated  on  the  tip 
vortex,  not  on  the  blade.  Thus  for  computing  the  distortion  of 
the  tip  vortex,  a  single,  rigid  wake  inboard  trailing  vortex, 
typically  trailed  from  the  50  percent  spanwise  station,  repre¬ 
sents  the  inboard  trailing  vortex  sheet  about  as  well  as  several 
rigid  wake  line  vortices  trailed  from  various  inboard  stations 
with  a  combined  circulation  equal  to  that  of  the  tip  vortex. 

To  compute  the  geometry  of  the  tip  vortex,  the  velocities 
induced  on  the  tip  vortex  by  itself,  and  by  the  inboard  vortex. 


must  be  calculated  at  various  times.  This  computation  involves 
the  integration  of  the  Biot-Savart  relation  over  the  tip  vortex 
and  the  inboard  vortex.  Since  this  integral  cannot  be  solved 
analytically  even  for  the  relatively  simple  rigid  wake  case, 
either  the  geometry  of  the  tip  vortex  must  be  approximated  by 
some  analytically  integrable  geometry  or  numerical  integration 
must  be  used.  As  shown  in  Ref.  7,  the  best  method  for  integrating 
the  Biot-Savart  relation  in  nonrigid  wake  cases  is  to  break  the 
trailing  vortex  lines  up  into  a  set  of  short  straight  vortex  line 
segments  each  subtending  an  equal  A<j>  (typically  A<{>  =  15  deg), 
where  <j>  is  the  azimuth  angle.  The  shed  vorticity  can  then  be 
represented  by  straight  vortex  line  segments  connecting  the 
corners  where  two  trailing  vortex  line  segments  meet,  both  on 
the  tip  vortex  and  the  inboard  vortex  at  the  same  azimuth  angle 
(<f>).  Figure  1  shows  two  turns  (4n)  of  the  spiral  rigid  wake  of 
one  blade  projected  into  the  tip  path  plane  (TPP),  where  £  =  1 
is  the  tip  vortex  and  £  =  0.5  is  the  inboard  vortex  trailed  from 
the  50  percent  radius  station. 

The  circulation  of  the  various  elements  of  vorticity  in  the 
model  must  be  defined  for  the  computation  of  induced  velocities. 
The  strength  of  the  bound  vortex  as  a  function  of  azimuth  angle 
(time)  is  determined  by  a  rotor  airloads  calculation  using  either 
a  distorted  wake  from  a  previous  iteration  of  the  wake  geometry 
computation  or  a  rigid  wake.  The  trailing  wake  vortex  line  seg¬ 
ment  between  azimuth  angles  <p  -  A<p  and  <p  is  taken  to  have  a  con¬ 
stant  circulation  equal  to  the  circulation  of  the  bound  vortex  at 
azimuth  angle  <f>.  The  shed  wake  vortex  line  segment  at  azimuth 
angle  $>  is  taken  to  have  a  constant  circulation  equal  to  the  dif¬ 
ference  between  the  circulation  of  the  trailing  wake  vortex  line 
segments  between  azimuth  angles  <}>  -  A<)>  and  and  between  azimuth 
angles  <t>  and  <j>  +  A4> . 

At  low  advance  ratios  (p  =  flight  speed/rotor  tip  speed) 
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when  the  effect  of  the  self-induced  velocities  on  the  wake  geome¬ 
try  is  greatest,  the  bound  circulation  has  a  large  steady  state 
value  compared  to  its  higher  harmonics.  A  first  approximation 
for  wake  geometry  calculations  (as  well  as  airload  calculations, 
see  Ref.  2)  is  to  assume  a  constant  value  of  the  bound  circulation 
based  on  the  mean  thrust  of  the  rotor.  This  results  in  a  trailing 
wake  where  all  of  the  vortex  line  segments  have  equal  circulation 
and  eliminates  the  shed  wake. 

This  vortex  line  representation  is  inadequate  when  it  is  ' 
necessary  to  calculate  the  induced  velocity  at  a  point  very  close 
to  or  on  the  tip  vortex.  For  this  case  the  effects  of  a  finite 
vortex  core,  where  viscosity  becomes  important  and  the  Biot- 
Savart  relation  no  longer  applies  directly,  must  be  included.  A 
useful  first  approximation  to  the  induced  velocity  inside  a  vor¬ 
tex  core  is  to  assume  solid  body  rotation  of  the  vortex  core  with 
an  induced  velocity  at  the  outer  edge  of  the  vortex  core  equal  to 
the  result  given  by  the  Biot-Savart  relation  at  that  point  (see 
Ref.  8).  Thus  for  most  cases  it  is  only  necessary  to  check 
whether  or  not  the  point  at  which  the  induced  velocity  must  be 
calculated  is  inside  the  vortex  core  and  apply  either  the  Biot- 
Savart  relation  or  the  solid  body  rotation  model. 

For  the  case  of  calculating  the  velocity  induced  at  a  point 
(P)  on  the  tip  vortex  by  the  two  immediately  adjacent  vortex  seg¬ 
ments,  the  use  of  straight  line  vortex  segments  is  inadequate, 
since  it  would  give  zero  induced  velocity  at  point  P  due  to  these 
two  vortex  segments.  These  two  vortex  segments  adjacent  to  point 
P  are  therefore  replaced  by  a  circular  arc  vortex  segment  passing 
through  the  three  points  determined  by  the  ends  of  the  vortex 
segments  (see  Fig.  A*0 .  The  velocity  induced  at  point  P  by  the 
circular  arc  vortex  segment  is  given  by  Ref.  9  as  a  function  of 
the  vortex  core  size,  the  radius  of  curvature,  the  circulation, 
and  the  vortex  segment  size  (A<J>). 
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Various  sources  estimate  the  vortex  core  diameter  at  from 
0.1  chord  to  1.0  chord.  Work  is  under  way  at  MIT  to  refine  the 
estimate  of  the  vortex-core  diameter,  and  to  define  the  v,oll-up 
process  in  more  detail.  For  the  case  of  a  rotor  in  forward 
flight,  considered  here,  the  results  turn  out  to  be  very  in¬ 
sensitive  to  vortex  core  size.  For  other  cases,  such  as  a 
highly  loaded  hovering  propeller,  the  details  of  the  rolling  up 
of  the  tip  vortex  just  after  it  is  trailed  are  likely  to  be  much 
more  important. 

Since  the  induced  velocity  calculated  by  the  Biot-Savart 
relation  falls  off  as  l/(distance)  vortex  elements  far  away 
from  either  the  rotor  or  the  portion  of  the  wake  being  distorted 
contribute  very  little  to  the  induced  velocity.  For  airload  cal¬ 
culations  only  the  first  two  or  three  turns  (2ir  in  azimuth  angle 
=  1  turn)  of  the  wake  are  significant  so  that  the  wake  geometry 
need  only  be  calculated  for  this  part  of  the  wake.  For  wake  dis¬ 
tortion  calculations  only  two  or  three  turns  of  the  wake  below 
the  wake  element  whose  distortion  is  being  calculated  plus  any 
wake  above  (nearer  the  rotor)  the  wake  element  need  to  be  con¬ 
sidered. 

Consider  the  vortex  element  at  the  lower  end( farthest  away 
from  the  rotor)  of  the  section  of  wake  whose  wake  geometry  is 
being  calculated  (two  or  three  turns  below  the  rotor).  To  cal¬ 
culate  the  distortion  of  this  vortex  element,  two  or  three  turns 
of  the  wake  below  this  vortex  element  must  be  considered.  For 
this  additional  section  of  the  wake,  an  extrapolated  version  of 
the  wake  geometry  must  be  used  since  its  wake  geometry  cannot  be 
calculated  directly.  The  extrapolation  procedure  used  here 
assumes  that  the  wake  has  the  same  distortion  parallel  to  the  tip 
path  plane  as  it  had  at  the  lower  limit  of  the  calculated  wake 
geometry  and  drifts  downward  perpendicular  to  the  tip  path  plane 
with  the  fully  developed,  uniform  rotor  inflow  velocity  from 
momentum  theory. 
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III.  DEVELOPMENT  OP  THE  COMPUTATIONAL  PROCEDURE 

Having  arrived  at  a  model  for  the  vortex  wake  of  the  rotor, 
the  calculation  of  the  wake  distortion  due  to  self-induced  ve¬ 
locities  can  be  considered.  A  straightforward  way  to  do  this  is 
to  start  with  M  (say  M  =  2  or  3)  turns  (2irM  in  azimuth)  of  rigid 
wake  and  compute  the  velocity  induced  by  this  wake  on  the  element 
of  the  tip  vortex  which  is  just  being  trailed  from  the  blade  tip. 
Advance  the  blade  in  azimuth,  generating  a  now  vortex  line 
segment,  and  compute  the  induced  velocity  at  both  the  point  where 
it  was  calculated  previously,  which  is  now  A<j>  behind  the  blade, 
and  at  the  tip  vortex  element  which  is  just  being  trailed  from 
the  new  blade  position.  Note  that  the  calculation  of  the  induced 
velocity  at  the  point  A<)>  behind  the  blade  only  requires  the  ad¬ 
dition  of  the  velocity  induced  by  the  newly  generated  vortex  seg¬ 
ment  to  that  previously  calculated  for  the  rest  of  the  wake.  The 
average  of  the  two  induced  velocities  calculated  at  the  point  A<f> 
behind  the  blade  can  be  multiplied  by  the  time  interval  A<|>/n  (where 
ft  is  the  rotor  angular  velocity)  to  obtain  the  wake  distortion  at 
this  point  during  this  time  interval. 

The  blade  is  again  advanced  by  a and  the  induced  velocity 
is  calculated  at  three  points:  at  the  blade  tip,  A<}>  behind  the 
blade,  and  2 A<f^  behind  the  blade.  Since  the  point  2 A«f>  behind  the 
blade  has  moved  relative  to  the  rest  of  the  wake,  the  induced 
velocity  at  this  point  must  be  calculated  by  integrating  over  the 
entire  M  turns  of  wake  again.  The  distortion  of  the  two  points 
A 4>  and  2A<J>  behind  the  rotor  is  calculated  as  before. 

This  process  of  advancing  the  blade  A$  and  calculating  the 
induced  velocity  and  hence  the  wake  distortion  at  an  ever  increas¬ 
ing  number  of  points  continues  until  the  distortion  is  being  cal¬ 
culated  over  m  <_  M  turns  (27rm  in  azimuth)  o*f  the  wake.  Then  as 
a  new  vortex  segment,  and  hence  point  at  which  wake  distortion  is 
calculated,  is  added  at  the  beginning  of  the  wake,  the  calculation 
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of  the  distortion  of  the  point  2-nm  behind  the  rotor  is  stopped 
and  the  extrapolation  procedure  is  applied  as  needed.  This  is 
because  the  wake  geometry  is  only  needed  for  the  first  m 
(typically  m  =  2)  turns  (27rm)  of  the  wake,  as  discussed  aoove. 

The  blade  continues  to  be  advanced  by  A<j>  and  the  wake  distortion 
of  the  first  m  turns  of  the  wake  calculated  until  a  stable  wake 
geometry  is  arrived  at. 

While  the  process  described  above  is  easily  programmed  to 
be  run  on  a  digital  computer,  it  requires  a  great  deal  of  compu¬ 
tation.  For  this  reason,  the  above  approach  has  not  been  tried  at 
MIT.  Instead,  the  initial  approach  assumed  that  the  wake  distor¬ 
tion  from  a  rigid  wake  was  small  and  an  iterative  procedure  was 
chosen.  This  starts  out  just  like  the  procedure  outlined  above 
but  when  the  wake  distortion  is  computed,  it  is  merely  stored  ana 
does  not  affect  later  induced  velocity  calculations.  For  the  en¬ 
tire  first  iteration,  a  rigid  wake  is  used  for  induced  velocity 
calculations.  When  the  blade  has  advanced  through  2um  in  azi- 
mut.i  ($),  the  first  iteration  ends.  The  wake  distortion  is 
calculated  from  the  time  history  of  induced  velocity  obtained  from 
the  first  iteration  and  is  used  to  calculate  a  wake  geometry  which 
is  then  used  for  all  induced  velocity  calculations  in  the  second 
iteration.  This  process  continues  until  convergence  is  achieved. 

This  iterative  procedure  saves  a  great  deal  of  computation 
over  the  "run  until  stable’1  method  since  the  points  at  which  in¬ 
duced  velocity  is  calculated  do  not  move  relative  to  the  rest  of 
the  wake.  Thus,  the  integration  of  the  Biot-Savart  relation  over 
the  entire  wake  does  not  have  to  be  done  repeatedly  for  the  same 
point.  Recall  that  "the  entire  wake"  means,  for  all  practical 
purposes,  only  M  turns  (2ttM)  of  the  wake  below  the  point  at  which 
induced  velocity  is  calculated  and  all  of  the  wake  above  that 
point,  as  discussed  previously. 

In  order  to  obtain  converged  solutions  with  this  technique, 
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some  refinements  are  needed.  This  is  because  of  the  [l/(distance )  ] 
variation  of  the  Biot-Savart  relation  for  induced  velocities.  When 
an  element  of  the  tip  vortex  finds  itself  in  a  region  of  high  in¬ 
duced  velocity,  it  is  swept  along  rapidly  with  the  induced  velocity 
and  is  likely  to  find  itself  soon  in  a  region  with  a  considerably 
different  induced  velocity.  In  a  method  which  only  allows  the 
wake  to  move  relative  to  itself  at  the  end  of  each  iteration  an 
element  of  vorticity  which  happens  to  be  in  a  region  of  high  in¬ 
duced  velocity  initially  will  remain  there  for  the  whole  iteration 
and  thus  accumulate  a  large  (and  unrealistic)  distortion. 

An  obvious  way  to  improve  this  situation  is  to  update  the 
wake  geometry  being  used  to  calculate  induced  velocities  with  the 
current  calculated  wake  geometry  every  few  (n^)  A<f> .  Thus  every 
tjme  the  blade  advances  through  n^A<j>  in  azimuth  the  Biot-Savart 
relation  must  be  integrated  over  the  entire  wake  again,  since 
the  points  at  which  induced  velocities  are  calculated  have  moved 
relative  to  the  rest  of  the  wake.  The  key  questions  in  this 
method  are  now  how  large  can  n1  be  and  how  many  iterations  are 
required  for  convergence  (see  the  Results  section  of  this  report). 
This  iterative  method  is  our  current  approach  for  calculating  wake 
geometry.  Appendix  A  outlines  the  computational  procedure  in  con¬ 
siderable  detail. 
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IV.  RESULTS 


Some  preliminary  results ,  run  to  evaluate  the  relative  im¬ 
portance  of  various  parameters  and  approximations,  are  presented 
here.  These  preliminary  results  were  all  run  for  the  same  stand¬ 
ard  case  (p  =  0.1125,  X  =  0.0*0  to  allow  them  to  be  compared  with 
some  experimental  wake  geometry  results  obtained  by  emitting 
smoke  from  the  tip  of  one  blade  of  a  two-bladed  rotor  running  in 
a  wind  tunnel.  The  wake  geometry  results  are  presented  in  the 
figures.  Airloads  have  been  computed  using  some  of  these  theo¬ 
retical  wake  geometry  results  but  they  are  not  presented  here  due 
to  a  lack  of  experimental  airload  data  for  comparison. 

The  experimental  (test)  data,  presented  here  is  from  some 
preliminary  work  to  be  reported  in  Ref.  11.  This  work  is  being 
continued  under  Army  Research  Office  sponsorship  (contract  DA-31- 
124-AR0-D-470) .  The  test  data  was  obtained  by  ejecting  hot  oil 
smoke  downstream  from  one  blade  tip  of  a  two-bladed,  articulated 
rotor  operating  in  the  return  section  (10  ft  x  12  ft)  of  a  closed 
circuit  wind  tunnel.  Simultaneous  photographs  from  below  and 
from  one  side  were  taken  and  the  coordinates  of  various  points 
on  the  tip  vortex  were  computed  from  this  data  with  the  aid  of  a 
high  speed  digital  computer.  These  coordinates  have  been  plotted 
and  are  the  curves  labeled  "test"  or.  the  figures.  Due  to  experi¬ 
mental  difficulties,  no  data  points  were  obtained  in  the  dashed 
portion  of  the  wake  and  this  has  been  sketched  in.  Although  the 
theoretical  wake  geometries  are  computed  for  a  full  4n  of  the 
wake  only  3*  to  3 .5"  of  test  data  have  been  obtained.  The  measured 
characteristics  and  operating  conditions  of  the  rotor  for  the 
standard  case  are:  R  =  2.425  ft.,  b  =  0.125  ft.,  o  =  O.O657,  fi  = 

200  rpm,  Locke  Number  =  2pubRVl  -  2.2,  V  =  353  fpm,  e(pitch  angle) 

=  7.1  deg,  i  =  -  12.5  deg,  a(shaft  angle)  =  15*1  deg.  The  cal¬ 
culated  rotor  operating  conditions  for  the  standard  case  are: 
v  =  0.1125,  CT  =  0.0036,  X  =  0.040,  -  u  tan  i  =  0.025,  yq  =  0.0465. 
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Due  at  least  partly  to  wind-tunnel  turbulence,  the  experi¬ 
mental  wake  geometry  has  not  been  found  to  be  completely  repeatable. 
Figure  2  shows  a  top  view  (in  the  TPP)  and  a  side  view  of  the  tip 
vortex  of  one  blade  of  a  two-bladed  rotor  for  two  different  ex¬ 
perimental  wake  geometries  and  rigid  wake,  all  for  the  standard 
case.  It  can  be  seen  that  the  tv/o  typical  test  results  have  the 
same  general  character  but  differ  in  detail.  One  of  these  test 
r  suits  is  used  for  comparison  with  the  theoretical  wake  geometry 
in  Figs.  3  through  8.  When  comparing  theory  and  experiment  it  will 
be  necessary  to  look  for  the  same  kind  of  agreement  in  general 
character  rather  than  complete,  detailed  agreement  because  of  this 
uncertainty  in  the  experimental  results. 

The  theoretical  calculations  whose  results  are  presented 
here  (Figs.  3  through  8)  have  all  used  m  =  2,  M  =  3,  and  A<J>  =  15 
degs,  chosen  from  experience  with  previous  wake  geometry  programs. 
These  results  are  intended  to  illustrate  the  effects  of  various 
approximations  and  parameters  on  the  calculated  wake  geometry. 
Figures  3  through  8  all  compare  a  theoretical  result  with  an  ex¬ 
perimental  (test)  result  and  with  a  X  =  0.0*1  rigid  wake. 

It  should  be  noted  that  while  the  figures  all  show  the  wake 
geometry  when  the  blade  is  at  ^  =  3**5  degs  for  comparison  with 
experiment, the  geometry  changes  as  the  blade  takes  different 
azimuth  positions.  Thus  a  complete  description  of  the  wake  geo¬ 
metry  would  require  a  separate  figure  for  every  blade  azimuth 
position  0  <_  ip  <_  2v .  The  computer  calculation  necessarily  gen¬ 
erates  the  information  needed  to  construct  such  plots  every  Atp 
(15  deg)  but  they  are  not  presented  here. 

Figures  3,  *1,  and  5  show  the  effect  of  varying  n1  on  the 
calculated  wake  geometry.  They  all  represent  essentially  con¬ 
verged  solutions,  use  a  nondimensional  vortex  core  radius  p  = 

0.05,  neglect  bound  vorticity,  and  assume  a  constant  strength 
trailing  wake  (yt  =  1.0),  which  implies  no  shed  wake  (yg  =  0.0). 
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Figure  3  presents  the  results  of  the  fourth  iteration,  n^  =  2j 
Fig.  4  presents  the  results  of  the  fifth  iteration,  n^  =  6,  and 
Fig.  5  presents  the  results  of  the  fifth  iteration  n-^  =  12.  Note 
that  the  n^  =  2  case  required  one  less  iteration  for  convergence. 
Convergence  is  determined  by  comparing  the  plotted  results  of 
successive  iterations  and  hence  is  somewhat  subjective.  The  total 
running  time  on  the  MIT  Computation  Center's  IBM  709*J  to  produce 
the  final,  converged  results  was  approximately: 


nl 

2 

6 

!  12 

Computer  Time  (min.) 

56 

25 

i  1* 

Comparing  Figs.  3,  and  5  it  is  seen  that  n^  =  2,  and  n^  =  6 
give  very  similar  results,  especially  in  the  side  view,  while 
n1  =  12  gives  somewhat  different  results,  notably  in  the  side 
view.  In  the  n1  =  12  case  part  of  the  wake  is  predicted  to  be 
much  closer  to  the  TPP  than  it  should  be  or  than  predicted  by 
n1  =  2  and  n.^  =  6.  Thus  n1  =  6  is  recommended  to  obtain  acceptable 
accuracy  for  most  purposes  while  requiring  less  than  1/2  as  much 
computer  running  time  as  the  n^  =  2  case. 

It  might  be  argued  that,  since  the  tip  vortex  is  by  far 
the  most  important  part  of  the  vortex  wake,  perhaps  the  rest  of 
the  vortex  wake  can  be  neglected,  at  least  for  the  purpose  of 
computing  the  geometry  of  the  tip  vortex.  This  has  been  tried  and 
Fig.  6  shows  the  result  of  the  third  iteration,  n,  =  2,  p  =  0.05, 

i.  V 

y.  =  1.0,  y  =  0.0,  with  bound  vorticity  neglected,  and  with 
i  -  1.0  only  instead  of  the  usual  l  -  1.0  and  0.5.  This  result  has 
not  converged  and  comparison  with  the  second  iteration  (not  pre¬ 
sented  here)  seems  to  indicate  that  it  is  diverging  at  least  over 
the  part  of  the  tip  vortex  farthest  from  the  blade  in  azimuth. 

In  any  case  the  results  are  sufficiently  different  from  the  i  = 

1.0  and  0.5  case  (Fig.  3)  to  indicate  that  the  inboard  (£.  =  0.5) 
vortex  cannot  be  neglected. 

Figure  7  shows  the  results  of  the  second  iteration,  n^  =  6, 
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p  =0.05,  including  bound  vorticity  and  with  the  y  and  y  de- 
termined  by  a  rigid  wake  airloads  calculation.  Note  that  conver¬ 
gence  only  required  two  iterations  as  compared  to  five  for  the 
Y  =  1.0,  y  =  0.0  case.  This  is  believed  to  be  due  to  the  damp- 
ing  effect  of  the  shed  wake.  The  total  running  time  to  conver¬ 
gence  on  an  IBM  709*1  was  approximately  13  minutes  plus  less  than 
one  minute  to  compute  the  y.  and  y  ♦  Thus  because  of  its  much 

U  S 

faster  convergence,  this  case  requires  less  computation  to  obtain 
a  converged  result  than  the  y+.  =  1*0,  yc  =  0.0  case  in  spite  of  . 
che  greater  computation  required  for  each  interation.  Comparison 
of  this  result  (Fig.  7)  with  the  y*.  =  1*0,  y  =  0  case  (Fig.  *1) 
shows  similar  agreement  with  the  test  results. 

Figure  8  presents  the  results  of  the  same  case  as  Fig.  7 

but  with  p  =  0.005  instead  of  p  =  0.05  to  check  the  effect  of 
c  c 

different  tip  vortex  core  radius  estimates  on  the  calculated  wake 
geometry.  Comparison  of  these  figures  shows  negligible  difference, 
which  is  fortunate  since  there  is  about  an  order  of  magnitude  un¬ 
certainty  in  p  . 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  work  reported  here  is  continuing  and  this  is  a  pre¬ 
liminary  report.  At  present,  a  reasonable  model  of  the  vortex 
wake  has  been  developed,  a  method  of  computing  the  geometry  of 
the  tip  vortex  has  been  developed,  and  the  relative  importance  of 
various  parameters  and  approximations  used  by  this  method  has 
been  evaluated  for  one  test  case  including  comparison  with  an  ex¬ 
perimentally  determined  tip  vortex  geometry.  As  a  result  of  this 
work  some  preliminary  conclusions  can  be  drawn  and  some  sugges¬ 
tions  for  further  refinement  of  the  vortex  wake  model  and  the 
wake  geometry  computation  procedure  can  be  made. 

1)  Although  the  tip  vortex  is  the  dominant  feature  of  the 
vortex  wake,  s;me  representation  of  the  inboard  trailing 
vorticity  must  be  included,  .nis  representation  of  the 
inboard  trailing  vorticity  can  be  fairly  crude  (for 
example:  a  single  vortex  line  trailing  from  the  50  per¬ 
cent  span  point)  and  if  it  is,  then  the  use  of  an  a  priori 
specification  of  its  geometry,  such  as  the  rigid  wake  ap¬ 
proximation  is  entirely  adequate.  A  possible  improve¬ 
ment  over  a  single  vortex  trailing  from  the  50  percent 
span  point,  which  has  not  been  tried  yet,  would  be  a 

set  of  rectangular  trailing  vortex  sheets  trailed  at 
constant  (spanwise)  strength  between,  for  example,  the 
30  percent  span  point  and  the  70  percent  span  point. 

The  downstream  dimension  of  the  rectangles  might  be 
taken  as  0.5R  A<j>  or  0.7R  At.  This  would  imply  some 
overlapping  and/or  gaps  between  the  vortex  sheets 
but  it  should  still  be  a  better  model  than  a  single 
trailing  vortex  line  without  requiring  as  much  extra 
computation  as  multiple  trailing  vortex  lines. 

2)  The  effects  of  a  time  varying  trailing  wake  circula¬ 
tion  and  the  corresponding  shed  wake  should  be  included. 
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The  current  model  which  accumulates  the  shed  wake  in  a 
series  of  line  vortices  joining  the  corners  of  the 
trailing  wake  model  (see  Pig.  1)  is  adequate  but  a 
possible  improvement , which  has  not  been  tried  yet, 
would  be  a  set  of  rectangular  shed  vortex  sheets 
similar  to  the  trailing  vortex  sheets  mentioned  above. 

3)  The  calculated  vortex  wake  geometry  is  very  insensi¬ 
tive  to  the  estimated  vortex  core  size  and  any  vortex 
core  radius  from  1.0  to  0.1  chord  gives  essentially 
the  same  results. 

4)  For  the  updating  interval  used  in  the  wake  distortion 
calculation,  n^  (see  Page  10),  a  value  of  n^  =  6  is  a 
good  compromise  between  computation  time  and  accuracy. 

A  great  deal  remains  to  be  done.  Although  the  current 
computer  program  requires  less  than  15  min  running  time  on  an 
IBM  7094  for  a  two-blade  case,  considerable  improvement  in  this 
respect  is  still  required  before  extensive  application  of  this 
technique  becomes  practical.  A  promising  way  of  improving  the 
current  computational  procedure,  which  is  now  being  tried,  is  to 
revise  the  method  of  determining  when  the  wake  distortion  used 
for  computing  induced  velocities  (D**)  is  updated  to  the  latest 
calculated  wake  distortion  fD).  This  is  currently  done  at  fixed 
times  whenever  the  %  blade  advances  n^A^  (see  Pages  10  and  36). 

It  is  proposed  to  keep  track  of  the  distortion  accumulated  at 
each  azimuth  angle  (<J>)  since  the  last  updating  of  D*  ’  for  that  \ J>. 
Whenever  the  magnitude  of  this  accumulated  distortion  exceeds 

some  limit  (D  )  then  D* 1  at  that  4>  only  would  be  updated  to  the 

_v  L 

latest  D.  D  is  a  parameter  similar  in  function,  but  not  in  magni- 
L 

tude,  to  n.^  in  the  current  program.  The  best  value  to  use  for 

D  will  have  to  be  determined  by  trial  in  a  manner  similar  to  the 
L 

determination  of  n^  =  6  as  a  good  compromise  in  the  Results  sec¬ 
tion  of  this  report.  Just  how  much  reduction  in  running  time  will 
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result  from  this  new  refinement  .‘mains  to  be  seen  but  it  is  ex¬ 
pected  to  be  considerable. 

An  alternate  way  to  reduce  the  computational  effort  re¬ 
quired  to  obtain  the  wake  geometry  and,  if  feasible,  the  most 
useful  approach  for  practical  applications  would  be  to  develop  a 
quasi-empirical  way  of  obtaining  an  a  priori  specification  of  the 
wake  geometry,  in  effect  a  better  rigid  wake  approximation.  This 
could  be  developed  by  building  up  experience  as  to  what  kind  of 
wake  geometries  can  be  expected  in  various  flight  conditions, 
using  both  experimental  data  from  flow  visualization  studies  and 
theoretical  results  obtained  by  computations  similar  to  those 
described  in  this  report.  Then  using  physical  reasoning,  plus 
perhaps  some  simple  analytical  approximations,  the  desired  speci¬ 
fication  of  the  wake  geometry  wou]  1  be  developed. 

Since  the  wake  geometry  is  only  important  to  the  extent 
that  it  influences  the  airloads  on  the  rotor  (or  fuselage)  a 
study  of  the  influence  of  the  wake  geometry  on  the  airloads  is 
necessary.  This  will  determine  just  how  good  a  quasi-empirical 
specification  of  the  wake  geometry  must  be  and  what  portions  of 
the  wake  geometry  are  the  most  important.  An  experimental  program 
where  the  airloads  on  a  rotor  are  measured  by  pressure  transducers 
at  the  same  time  as  the  geometry  of  the  tip  vortex  is  determined 
by  flow  visualization  would  be  very  useful  in  support  of  this 
effort . 

A  question  which  has  not  been  considered  here  but  which 
is  vital  to  a  full  understanding  of  rotor  airloads  and  noise  is  that 
of  shore  range  blade-vortex  Interaction.  Due  tc  wake  distortion, 
a  rotor  blade  occasionally  passes  very  near  to  or  intersects  the 
tip  vortex.  This  results  in  large  airload  peaks  and  considerable 
amounts  of  very  objectionable  noise.  It  may  also  cause  the 
breakdown  of  the  tip  vortex  locally  with  a  resultant  effect  on 
the  airloads  and  the  wake  geometry.  This  problem  is  exceedingly 
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complex  and  will  require  the  representation  of  the  blade  by  a 
lifting  surface  and  probably  the  consideration  of  both  viscosity 
and  Mach  number  effects.  The  ultimate  aim  of  work  on  this  prob¬ 
lem  should  be  a  good  enough  understanding  of  the  problem  to  allow 
the  formulation  of  simple  quasi-empirical  approximations  for  the 
airloads  and  noise  generated  by  these  interactions  and  to  indicate 
what  happens  to  the  vortex  after  such  an  interaction. 
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FIGURE  I:  MODEL  OF  THE  FIRST  4  v  OF  THE  WAKE  OF  ONE  BLADE  PROJECTED  INTO 
THE  TPP  (RIGID  WAKE) 


VIEW  IN  TPP 


FIGURE  2:  COMPARISON  OF  TEST  (FLOW  VISUALIZATION)  RESULTS 
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SIDE  VIEW 


- TEST  — THEORY  -  RIGID  WAKE 


FIGURE  4  COMPARISON  OF  THEORY  AND  TEST  RESULTS - 
n,=  6,5th  ITERATION,  pc  =  0.05,  /t  =1.0,  /S  =  0.Q 
AND  4=1.0, 0.5 
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FIGURE  6  COMPARISON  OF  THEORY  AND  TEST  RESULTS 
n ,  =  2 , 3 rd  ITERATION  ,  pc  =  0.05  ,  =  1.0  , 

/s  =  0.0  ,  AND  1=  1.0 


SIDE  VIEW 


TEST  — • —  THEORY  -  RIGID  WAKE 


FIGURE  7  COMPARISON  OF  THEORY  AND  TEST  RESULTS - 

n,  =  6,  2nd  ITERATION  ,/>c  =  0.05  ,  xt  AND  FROM 
RIGID  WAKE  AIRLOADS,  AND  l  =  1.0 , 0.5 


ROTOR  BLADE 


VIEW  IN  TPP 


TEST  —♦—THEORY  - RIGID  WAKE 


FIGURE  8  COMPARISON  OF  THEORY  AND  TEST  RESULTS  - 
n)  =6, 2nd  ITERATION  ,/>c=  0.005,  yx  AND  /«.  FROM 
RIGID  WAKE  AIRLOADS  ,  AND  2  =  1.0 , 0.5 


APPENDIX 


THE  COMPUTATIONAL  PROCEDURE 

The  vortex  wake  geometry  can  be  broken  up  into  two  parts: 
the  basic  geometry  due  to  the  motion,  relative  to  the  free  stream, 
of  those  portions  of  the  rotor  blades  from  which  the  wake  is 
trailed,  and  the  wake  distortion  due  to  the  velocity  field  in¬ 
duced  by  the  wake  and  the  bound  vorticity.  This  division  of  the 
wake  geometry  is  convenient  since  the  basic  geometry  is  easily 
calculated  for  a  given  advance  ratio  (u  =  V/flR)  and  tip  path 
plane  incidence  relative  to  the  free  stream  (i),  see  Fig.  Al, 
allowing  effort  to  be  concentrated  on  the  calculation  of  the  wake 
distortion. 

For  the  computation  of  the  nondimensional  wake  distortion 
(D)  of  point  P^  on  the  vortex  wake,  a  coordinate  system  is  chosen 
with  its  origin  at  the  point  where  P^  would  be  if  only  the  basic 
geometry  was  considered.  A  right-handed  axis  system  is  defined 
with  the  y  axis  and  j  parallel  to  and  in  the  opposite  direction 
from  the  line  of  rotor  hub  motion  and  parallel  to  the  tip  path 
plane  (TPP).  The  x  axis  and  i  are  perpendicular  to  the  y  axis  and 
are  parallel  to  the  TPP.  The  z  axis  and  k  are  perpendicular  to  the 
TPP  and  point  downward  (see  Fig.  Al).  i,  j,  and  k  are  unit  vectors. 

Let  the  time  at  which  the  wake  distortion  of  point  P^  is 

calculated  be  t^.  Let  the  time  at  which  the  element  of  vorticity 

at  point  P  was  trailed  be  t  .  Then  at  time  t  point  P  was 
n  n  n  n 

located  on  a  radial  spanwise  line  from  the  center  of  the  rotor 
hub  (representing  the  quarter-chord  line  of  the  blade),  at  a  dis¬ 
tance  nR  (n  =  1.0  for  the  tip  vortex),  from  the  hub  on  the  blade 
with  azimuth  angle  tj>(the  n  blade),  where  R  is  the  rotor  radius. 

Consider  a  typical  vortex  segment,  whose  induced  velocity 

at  P  must  be  calculated.  Let  the  time  at  which  the  last  element 
n 
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of  vorticity  forming  part  of  the  vortex  segment  was  trailed  be 
t  .  If  the  location  of  this  element  of  vorticity  is  labeled  P 
then  at  time  t^  point  P  was  located  on  a  radial  spanwise  line 
from  the  center  of  the  rotor  hub  (representing  the  quarter-chord 
line  of  the  blade),  at  a  distance  £R  from  the  hub  on  the  blade 
with  azimuth  angle  <j>(the  £  blade).  Note  that  the  first  element 
of  this  vortex  segment  was  trailed  at  time  t  —  A<|>/fl  when  the 
£  blade  was  at  azimuth  angle  (A  =  rotor  angular  velocity, 

radians/sec) . 

Let  the  instantaneous  angle  between  the  n  blade  and  the  £ 

blade  be  x,.  Then  t  —  t  =  +  ?  — <f> )/fl,  for  example. 

n  ~ 

The  distortion,  at  time  t^,  of  the  element  of  the  tip  vor¬ 
tex  located  at  point  P^  is  a  function  of  two  variables.  One  is 
the  azimuth  angle  (ijO  of  point  P^  at  the  instant  of  trailing  (t^) 
The  other  is  the  time  since  the  element  of  vorticity  was  trailed 
td  -  t  =  (<j>'  —  tfi  -  O/n,  where  <t>  •  is  the  azimuth  angle  of  the  £ 
blade  at  time  td»  The  range  of  (<J>*  —  <l>  -  c )  is  determined  by  the 
number  of  turns  of  the  wake  (m)  whose  distortion  is  calculated 
-  <!>-?)<_  2irm,  while  the  range  of  is  0  <_  <J>  <_  2tt. 

The  nondimensional  total  wake  distortion  ("D)  is  the  sum  of 
the  nondimensional  distortions  due  to  the  trailing  wake  (D  ), 

_ V  ^ 

the  shed  wake  (Ds),  and  the  bound  vorticity  (DgV). 

o=  +  + 

The  wake  distortion  due  to  the  trailing  wake  is  the  distor 
tion  which  accumulates  in  time  A$/n  due  to  the  vortex  trailed 
from  the  £th  spanwise  station  on  the  blade  summed  over  time, 
v'oth  £ ,  and  all  c . 
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The  summation  over  4>'*  (time)  represents  the  integration  of  the 

induced  velocity  over  time  while  the  n  blade  moves  from  azimuth 

angle  ip  at  time  t  to  azimuth  angle  <p '  —  c  at  time  t,.  This  is 

n  fl 

the  same  as  the  *  blade  moving  from  f  +  t,  to  <j>*. 

Similarly,  the  wake  distortion  due  to  the  shed  wake  is  the 
distortion  which  accumulates  during  time  A<{)/n  due  to  the  shed 

t* 

wake  trailed  by  the  5  blade  summed  over  time  and  all  c. 

Z  ADs(s,*(f) 

The  nondimensional  wake  distortion  due  to  the  bound  vorticity  is: 

0  t'-V-s)  =  Z  Z  ADt 

0V  <•  in  .  , 

Note  that  no  sum  over  z  is  needed  for  the  shed  wake  or  the  bound 

vorticity  since  for  the  21  case  considered  here  there  is  only  one 

element  of  shed  wake  or  bound  vorticity  at  any  azimuth  (see  Fig.  1). 

In  all  three  cases  above, the  sum  on  <p '  '  starts  at  ip  +  c  +  A<j> 

because  at  (j*1'  =  ip  +  ;  no  time  has  elapsed  since  P  was  trailed 

n 

and  hence  no  distortion  has  occurred.  Thus,  when  <J>'  =  +  c,  set 

D(ij»,  <j> '  —  ip  -  c)  =  0  for  all  ip(0  ^  £  2tt). 

Let  the  nondimensional  induced  velocity  at  point  at  the 
instant  when  the  Z  blade  has  azimuth  angle  <t> '  ’  due  to  the  vortex 
trailed  from  the  £th  spanwise  station  on  the  blade  be 
"q  ( z ,  5  ,  ip  ,<j> '  ' ) .  Assuming  that  the  average  Q  in  a  time  interval 
is  the  average  of  the  Q  at  the  beginning  and  end  of  that 
time  interval: 
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a\u,  s,  t  ?)  -kxk  [Oti,s,  t<n+ou,\ t  f- 

_  Q 

where  u  is  the  blade  semichord,  Ya~  2.TTy\fi.R  and  r0  is  the 
time  average  bound  circulation.  The  nondimensional  induced 
velocity  due  to  the  individual  vortex  segments  (cf^ )  is  summed  to 
obtain 

$ 

Q«, s,%*>  E  Y,u*>  itlJAt*) 

<f>=</'+S-2trM  +  A(t> 

where  W>“  rt(i»/C  and  rt  is  the  circulation  of  an 
individual  trailing  vortex  segment. 

Averaging  is  not  used  for  the  shed  wake,  therefore: 

f'.  A <j> 

t=t+S-Z-rM 


where  X  (4)  •  m/Q  ,  rs  is  the  circulation  and  cfg  is 
the  nondimensional  induced  velocity  of  an  individual  shed  vortex 
segment . 

For  the  bound  vorticity  no  sum  over  <j>  is  needed. 


From  Ref.  10,  the  induced  velocity  at  point  due  a  vor¬ 
tex  segment  of  circulation  r  is  perpendicular  to  the  plane  formed 
by  the  vortex  segment  and  point  and  has  the  magnitude 
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^  r  Vc°s  **+ 

,4r R  /\  fi/R  /  (see  Fig.  A2). 
nondimensionalized  by  (r/^uR)  is  then: 


The  induced  velocity 


(-a)-(-g)  ,  (-t)-c 

1*1  ID  / 

|(-c)x(-S)| 


see  Fig.  A2. 

Let  a,.,  b  .  and  c,  be  used  to  define  q  and  a.  b  .  and  c  to 

bSSS 

define  qs.  For  the  trailing  wake: 

xr-x(J,t)  ,  y  =  yUi  <t>)  t 

2=  x,  =  x(*,^-a4)  1  , 

and  p  ss  Jl  1 1  <i -£<&}.  Then  from  Fig.  A2: 

ct=  tt  -1t  .  For  the  shod  wake:  X=  X (1.0,4), 

a- y(i.o,4)  ,z.z(i.o,*)  x^xfo.s.4)  3,=  u^-s,*) 

a,  =  *(0.5,4).  '  '  1  ‘ 

Prom  Fig.  A2 :  ^  (,.0, 4)  *  +  +  a  (1.0,4)  W  , 

tj*  xfo.S, 4»> t  +3(0.5, 4)J  +  z(o.^,4)k  ,  cs=as-ts 

Let  the  nondimensional  distoition  of  point  be 

02=  DIZt+O^Zf+DKZk3^  that  of  P°lnt  pn  be 

53=  0I3t +  033J +  0K31J  - 

From  Fig.  A3: 

x(j,4)  =  -  *<"  ^  +  012  -  0X3  , 

3 (i,4)  =  i«s  4,  -  eostJ'+^^+s-^  +  o^-oaB  . 
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Recalling  that  the  tilt  of  the  tip  path  plane  relative  to  the 
free  stream  (Fig.  Al)  causes  a  displacement  of  the  wake  in  the 
z  direction  of  —  ptani  (  ^  s-4  )  similar  to  the  -<f>) 

term  in  the  y equation  (recall:  f  - 

zU,4')=-/i-twi(f+S-<l‘)  +  ^U-I.o)  +  DK2-OK3) 


where  aQ  is  the  coning  angle  (see  Fig.  Al). 


When 


;ne  coning  angle  (.s 

_h.s 

r  ici  A 


c  point  is  inside  the  vortex 


core,  where  p  R  is  the  radius  of  the  core  of  the  tip  vortex  and 

V  ^ 

a  modified  definition  of  q  is  needed  (see  Page  6  ).  First  the 

induced  velcity  (q  )  at  h/R  =  p  ,  the  edge  of  the  vortex  core 

p  c  c 

is  determined. 


h 


* — 

- 1 

i 


v<^  ' 


i 


ax  b 


£x  b| 


where 


and 


c,  **  |a|  cos  <  - 
c  =  I  b  I  /3  ~ 


ic| 

(rb)-- 

Icl 


The  Induced  velocity  (q)  at  a  point  inside  the  core  can  then  be 
written,  assuming  solid  body  rotation  of  the  vortex  core: 


f 


K-a)x(-*)| 
1*1  />« 


As  mentioned  in  the  Vortex  Wake  Model  section,  the  trailing 
vortex  segments  A<J>  on  either  side  of  point  P^  are  considered 
separately  from  the  rest  of  the  wake  to  include  accurately  the 
effects  of  vortex  core  size.  Reference  9  gives  an  expression 
for  the  induced  velocity,  due  to  a  circular  arc  vortex  segment  with 


a  finite  core  radius  (pcR)  and  pc  <<  1,  at  one  end  of  the  vortex 
arc.  Put  into  our  notation,  this  gives  an  expression  for 
when  5 -0i  X-l.O  ,  and  or 


where  ?,  =  tf  (l.O ,  +  C,-  Cf  (l.O ,  S)  } 


As  shown  in  Fig.  A4,  the  radius  of  curvature  (p)  used  is  that 
for  a  circular  arc  passed  through  the  three  end  points  (one  of 
which  is  common)  of  c^  and 

When  c  =  0  and  $  =  ip  +  c  point  is  located  on  a  shed 
or  bound  vortex  segment  and  since  a  straight  line  vortex  in¬ 
duces  no  velocity  on  itself  =  0. 

The  key  approximations  used  here  to  obtain  more  rapid  com¬ 
puter  solutions  are  those  involved  in  the  specification  of  D2  and 

D3  the  distortion  of  points  P.  and  P  ,  respectively.  When  the 

x*  n 

rigid  wake  approximation  is  used,  this  is  straightforward.  At 
time  td  point  P£  is  located, using  the  rigid  wake  approximation: 


is  the  nondimensional,  uniform  inflow  velocity,  perpendicular  to 
the  TPP,  from  momentum  theory  and  C,p  is  the  thrust  coefficient 
of  the  rotor.  At  time  t^  point  P^  is  located  using  the  rigid 
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wake  approximation: 

0X3  =  033=0,  DK3  = 

It  may  be  helpful  her 


to  recall  that: 


5) 


tjL-tx=  (</>'-  4)/ji 

and 

tjL~ty=  fol  tf'  -s)/j2. 

The  method  of  determining  D2  and  D3  for  the  general  case 
when  a  nonrigid  wake  approximation  may  be  used  is  best  illus¬ 
trated  by  an  outline  of  the  important  steps  in  the  wake  distor¬ 
tion  computation. 


1) 


Input  the  data  including  D‘,  which  is  the  distortion 
(D)  output  by  the  previous  iteration.  For  the  first 
iteration,  the  rigid  wake  approximation  is  used  to 
define  D' : 


01  =  03=0,  dk'= 


jut 


where  At  is  the  time  elapsed  since  the  element  of  vor- 

•najb 

ticity  whose  distortion  is  to  be  approximated  by  D' 
was  trailed  from  the  rotor  blade. 


2)  Defjne  D1'  as  the  distortion  used  to  determine  the 

wake  geometry  used  for  induced  velocity  calculations. 
Initially  set  D* '  =  D'. 


3)  Set  a  counter  n  =  1 

D)  Compute  0  (¥*,  S)  for  °  ~  t  ~  Zt 

and  nn.A* 

During  this  calculation  D'1  remains  unchanged  which 
implies  that  f’or  any  given  ^ 

combination  will  remain  unchanged  for  all  values  of 

*'  suoh  that  nv,A<f>  ■ 
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and : 


5) 

6) 
7) 


This  means  that: 

SU, *,*'**)+  f(*i*,*i*') 

5(l/',^,-i/'-s)=  <j>‘) 

Which  avoids  recomputing  all  of  the  ^  (Ji}% 
and  summing  to  get  and 

for  each  different  <j> ' . 

Now  update  D**  ' .  Let  *D"  j1-  #  -  S  )  =  0  f'-fl-S) 

for  O^^ZTTand  (n-l)n,  6  </>-  5)  ^  T1  ttj  A<f> 

If  T1  TJj  A<|>  > WI  end  of'  iteration,  otherwise  continued. 

n  =  n  +  1  go  to  step  4.  Now  because  D’ 1  was  changed  in 
step  5,  all  of  the  f  U,*M)  will  have  to  be  recomputed 
in  step  4 .  However,  they  need  only  be  computed  once 
for  each  time  step  4  is  reached  (see  the  explanation  of 
this  under  step  4).  From  this  description,  it  can  be 
seen  that  n.^, which  determines  how  often  D1’  is  updated, 
should  be  as  large  as  possible  (without  compromising 
convergence)  to  minimize  the  number  of  times  q*  must  be 
recomputed. 

—A 

Now  D2  and  D3  can  be  computed.  For  locating  point  P£ 
when  1  =  0.5j  the  rigid  wake  approximation  is  always  used. 

For  locating  point  P£  when  £.  =  1.0  and  for  locating  point  P^ 
at  all  times  D2  and  D3  are  calculated  in  terms  of  the  D1  and  Tr  1 
defined  above.  For  ,  D1  *  is  equal  to  *D  and 

should  be  a  better  approximation  to  the  distortion  than  D*  from 
the  previous  iteration,  therefore: 

D2*  0"  (+,♦-*)  for  £  (n-fln,  a4  . 

For  (n-l)Yl,  -  4  “  ^  ~  ^Tm  s  5*"  is  still  equal  to  D*’ . 
A  better  approximation  to  the  distortion  in  this  range  than  just 
using  D’  is  to  use  L’1'  for  the  distortion  accumulated  in  the 
initial  time  interval  and  add  the  D'  distortion 

Jl 
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accumulated  in  the  last 


time  interval; 


JT, 

therefore: 

02  *  ID 1 “  (<f>}  (n-On,  a<£) +  0  (<t ,  ♦*-  ♦)  -  ^  (' ♦  »^*1^  *i 

for  (Yi-i)n,A<j>  4  - 

For  2u*Vi1  4  the  distortion  is  not  calculated.  The 

distortion  for  this  range  is  approximated  by  using  the  definition 
above  for  the  distortion  accumulated  in  the  initial  ZlTYH^Sl.  time 
interval  and  the  rigid  wake  approximation  for  the  distortion 
thereafter.  Therefore, 

0X2  6$ ,$'-♦)  =  0X2  (<t>,2'ir»'i) 

03  2  (<M'-<fr)  a  032 

0K2  (4, 4>'-  4»)  =  DK2.  (4>,  2irm) 

for  2irm  <  4>'-  4>  ■  ^ 

since  (n-iln.A^  ift'-f-s)*  ,  03 

is  deigned  in  a  fashion  similar  to  D2  in  this  range: 

03=  D"(4>-ikA<t>)+  o'(^  f-t-s)-  o'0;Mv*). 

A  few  auxiliary  ca  oulations  are  needed  to  find  some  of 
the  inputs  for  the  computational  procedure  outlined  above.  The 
thrust  coefficient  (C,p)  is  defined  conventionally , 

-  X 

Ct”  fir f cm)1 
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where  the  thrust  ,  p  is  the  atmospheric  density, 

W  is  the  weight  of  the  aircraft  and  D  is  the  drag.  The  tip  path 
plane  incidence  (i)  can  be  estimated  by:  tan  i  *(d/w)  especially 
for  low  p  where  wake  distortion  is  most  important.  The  coning 
angle  (aQ)  can  be  estimated  from  uniform  inflow  theory  or  from 
a  rigid  wake  airload  calculation. 

The  time  average,  nondimensional ,  bound  circulation  (y  ) 
is  found  by  requiring  that  the  proper  thrust  (or  CT)  is  developed. 
The  lift  per  unit  span  is  L  =  pVT.  The  thrust  is 

where  n^  =  number  of  blades.  In  terms  of  time  average  values: 

U7>VoC  and  T=  JnfcS  L. 

'  o 

where 

ana  [>  2-irbJiRY. 

Recalling  the  definitions  of  and  o,  and  the  approximation 
that  there  is  constant  circulation  between  n  =  (normally 
0.5)  and  1.0  and  zero  circulation  inboard: 
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A  first  estimate  of  the  circulation  of  the  vortex  wake 


neglects  the  higher  harmonics  and  takes  the  time  average 
(zeroth  harmonic)  bound  circulation  (r  )  as  the  circulation  of 
the  trailing  wake.  This  implies  yt  =  1.0  and,  conserving  circu¬ 
lation,  y  =0.0. 

s 

An  improved  estimate  of  the  vortex  wake  circulation  can 
be  obtained  from  a  time  history  of  angle  of  attack  (a)  as  cal¬ 
culated  by  an  airload  calculation  (either  uniform  inflow,  rigid 
wake,  or  the  wake  geometry  from  a  previous  iteration  may  be  used 
in  this  airload  calculation).  The  angle  of  attack  ct(iji,n)  is 
computed  both  as  a  function  of  azimuth  angle  ^(tirne)  and  radial 
station  n.  The  lift  per  unit  span  is  L  =  pVr  or  alternately, 

L =  £oV*2bCu  and  CL  =  =  *  a.  . 

Since  V=J18  (vj  +  ^sinl/')  and  P  =  Zirb  JlR  Y 

L=  fljLR  (y+y*  *”> •#’)  Ztt  bjlR  Y 

L*>  />«(>  (J1R)Z  tyu  s!n</')Z^  . 


The  thrust  of  a  given  blade  as  a  function  of  azimuth  is 

A 


where  nQ  is  the  nondimensional  root  cutout. 


LLUJy  , 

Requiring  the  T^yO  developed  by  the  model  used  here,  which 
assumes  y(t^)  is  constant  from  n  =  Jtj  (normally  0.5)  to  n  =  1.0 
(the  blade  tip)  and  zero  inboard  be  equal  to  the  Tb(i|>)  computed 
by  the  airloads  calculation:  ^ 

ni  1  ‘  ‘  ZT 


'  ’  * 

/  (7+^ si" 


40 


Recalling  that  y,  is  normalized  by  y  : 

t  o 


» 

i 

ffyyusini *  O', y) 

h 

;i(i-jtx‘)+ ft--*,) 

Then 


and  conserving  circulation } 


X  (♦■)  -  vtd.o,  ♦+ Af)  -  y  . 


POINT 


FIGURE  Al:  TIP  VORTEX  BASIC  GEOMETRY  (  £  =0  CASE) 


FIG.  A2  VORTEX  SEGMENT  GEOMETRY 
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